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Generalized uncertainty relations are known to provide a minimal length h^/]3. The effect 
of such minimal length in the Casimir-Polder interactions between neutral atoms ( molecules ) 
is studied. The first order correction term in the minimal uncertainty parameter is derived 
and found to describe an attractive potential scaling as r~ 9 as opposed to the well known r~ 7 
(-} • long range retarded potential. 
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INTRODUCTION 

The complete understanding of the properties of the quantum vacuum is a point of central 



importance to both fundamental and applied physics [1]. A celebrated mechanical manifestation 



£SJ . of the fluctuations of the electromagnetic field around its zero point (vacuum) configuration is 

> ' Q 

the well known Casimir effect |2|] whereby an attractive force arises between two metallic neutral 

o : 

plates separated by a distance d . The Casimir effect has recently received a lot of attention mainly 



because of the fo 
measurements 



lowing reasons: i) advances in the experimental techniques leading to precision 



and to the possibility of measuring the dynamic Casimir effect [6j; ii) increased 
importance of the Casimir forces in the field of micro-electro-mechanical systems (MEMS) ?J and 
nano-devices [8(|; Hi) relevance to physics beyond the standard model of particle physics, such as, 
for example, hypothetical extra-dimensional models The present work is a contribution to the 
last of the above points and tries to make a connection between the Casimir effect and models 
which are generally refereed to as minimal length theories. 

The concept of a minimal length arises naturally in quantum gravity. Indeed, when trying to 
resolve small distances, higher energies are needed which eventually will affect the structure of 
space-time via their gravitational effects. Clearly this type of effect should occur at an energy scale 
of the order of the Planck mass. While it is generally stated that the incorporation of gravity within 
quantum field theory spoils the renormalizability of the latter, the fact that quantum gravity could 



1 See however ref. 0] for a different viewpoint which claims that there is no need to invoke vacuum fluctuations to 
explain the Casimir effect. 
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play an important role in the suppression of infinities in quantum field theories has been pointed 
out long ago [101 ] . 



Other models have been discussed in the literature which are also related to a minimal length 
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as 
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models might turn out to be accessible experimentally in the next generation of high-energy accel- 
erator such as the Cern LHC or the International Linear Collider. 

Although not trivial to prove it is natural to expect that a minimal distance should correspond 
at the quantum level to a minimal uncertainty. On the other hand string theories have been 
shown to predict the existence of a minimal length in the form of an effective minimal uncertainty 
(Ax) m i n |l4j |. This has triggered the study of quantum mechanical models based on generalized 
commutation relations such as: 

[x,p] = ih(i + pf) , (l) 

which lead to a generalized uncertainty principle (GUP) which in turn provides a minimal un- 
certainty. In particular in Ref. an Hilbert space representation of the generalized uncertainty 
relations is constructed explicitly in a one dimensional model, as well as in higher dimensions, iden- 
tifying the states which realize the maximal localization (i.e. the minimal uncertainty) in position 
space. An important aspect of these theories is that the eigen-states of the position operator are no 
longer physical states and so one is forced to introduce the so-called quasi-position representation, 
which consists in projecting the states onto the set of maximally localized states, instead of using 
the standard position representation obtained projecting the state vectors on the eigen-states of 
the position operator. Thus momentum states which are normally represented by plane waves 
(27r) _3//2 exp[(i/h)p ■ x] = (x\p) are now replaced by the set of functions tfjp IL (x) = (ip^ IL \p), 
where ip^ L are the maximally localized states (minimal uncertainty) around an average position 
x. The particular form of these states depends on the number of dimensions and on the specific 
model considered. These type of models are also connected to non-commutative quantum field 
theories. 

In Ref. [1(J the author studies the standard Casimir effect of two perfectly conducting plates at 
a distance d computing the correction term arising within a model based on generalized uncertainty 
relations: 
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It is worth to be noted that the correction term due to the minimal length (Ai) m j n ~ h\f]3 scales 
with a different power law (<i~ 6 ) and it is repulsive. This result is used to obtain an upper bound 
on the minimal length which turns out to be (Ax) m j n < 150 nm. 

In this work we take up the above scenario of generalized uncertainty relations and compute the 
correction term to the intermolecular Casimir-Polder interactions induced by a minimal length, 
expecting to find results similar to those of Ref. [lf| (c.f. Eq. [2]). In particular, by performing 
a quantization procedure of the electro-magnetic field as in Ref. [lfl, it will be shown that if a 
minimal length exists in nature the retarded Casimir-Polder interactions between neutral atoms 
and molecules will acquire, in addition to the standard r~ 7 interaction, a corrective term which 
scales as r~ 9 . However, as opposed to the plate-plate case, c.f. Eq. ([2]), the new term has the same 
sign of the standard QED result, i.e. it describes an attractive interaction. 

The reminder of the paper is organized as follows. In Section [TT] we discuss the generalized 
uncertainty relations and define a set of maximally localized states; in section IIIII we discuss the 
quantization of the electromagnetic field in the presence of a minimal length and, in section ITVl we 
derive the corrections to the Casimir-Polder intermolecular interactions due to a minimal length. 
Finally in section [V] we present the conclusions. 



II. GENERALIZED UNCERTAINTY PRINCIPLE 



Let us consider the generalized commutation relations of Eq. [TJ for simplicity in one dimen- 
sion. From these modified commutation relations one derives the generalized uncertainty principle 

(GUP) Q 

AxAp>^[l+P(Ap) 2 + (p} 2 ] p>0, (3) 

which is found to be related to a minimal length [if], [rJ. Indeed with a simple minimization 
procedure it is easily verified [15] that Eq. ([3]) implies an absolute minimal uncertainty (Ax) m j n = 
Hy/]3, which is obtained for those states such that (p) = and (Ap) 2 = (p 2 ) = 1/(3. An Heisenberg 
algebra that satisfies the generalized commutation relations of Eq. (JTJ) is represented on momentum 
space wave functions by: 

ptp(p) = p^{p), (4) 
x ip(p) = ih(l + (3p 2 ) d p ip(p) . (5) 

In Ref. [ljj the eigenstates of the position operator have been shown to be non physical because 
their uncertainty in the position (which vanishes) is smaller than the absolute minimal uncertainty 
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(or minimal length). One then can look for a set of maximally localized states for which the 



uncertainty in the position is minimal (i.e. equal to the minimal length) [17]. The procedure 
proposed by Kempf, Mangano and Mann (KMM) 15] to find the maximally localized states consists 
in minimizing the value (Ax)^ between those states which realize the generalized uncertainty 
principle in Eq. ([3]) with the equality sign (i.e. a squeezed state): 



(Ax) i ,{Ap) i , = - | ([x,p])^\ 



(6) 



In the adopted momentum space representation Eq. 



'mi + m d p - ( X ) + lh i±E^l±m {p _ m 



takes the form of a differential equation: 

lfr(p) = 0, (7) 



which admits the normalized solution ((a 



TT 



2(Ap) 2 
e, (p)=0, Ap=l/VP): 

£ arctan(y^p) 



(1 + Pp 2 ) 2 exp 



(8) 



As opposed to ordinary quantum mechanics it turns out that such states, in addition to being 
normalizable, are of finite energy and no longer "orthogonal" i.e. their closure relation involves a 
finite function instead of a Dirac-<5 distribution — £')): 
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Extension to higher dimensions. In Ref.[l7] the above ideas have been extended to a number n 
of arbitrary (spatial) dimensions. An important fact to recall is that there is no unique extension 
of Eq. [1] in more than one dimension. Indeed in order to preserve translational and rotational 
invariance generalized commutation relations must take the form 

[xi,pj] = ih [f{p 2 )dij + g(p 2 )piPj] , ij = l,...,n. (10) 

The functions f(p 2 ) and g(p 2 ) are not completely arbitrary. Relations between them can be found 
by imposing translational and rotational invariance. In the following we take up the standard choice 
of Ref. [itJ which has become rather popular and has been the object of many phenomenological 



studies: 



f{f) = —J^= , g{f) = (3. (11) 
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A spectral representation can then be found such that: 



(12) 



The KMM procedure consists in minimizing the position uncertainty within the set of squeezed 
states, see Ref. P, Q for details, and yields the following maximally localized states around a 



mean position £ (in the spectral representation): 



ip¥ L (z) = N [1 — P^— 



2\ Q / 2 



exp 
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(13) 



where a = 1 + v/T+ n/2 ~ 2.58 (in three spatial dimensions) is a numerical constant that charac- 
terizes the maximally localized states in the KMM approach. On the other hand when the number 
of dimensions n > 1 the generalized uncertainty relations are not unique and different models (ac- 
tually an infinite number of them) may be implemented [17] by choosing different functions f(p 2 ) 
and/or g(p 2 ) (c.f. Eq. [TT]) which will yield in general different maximally localized states which 
will not contain at all the numerical constant a. The minimal position uncertainty for the KMM 
maximally localized states is 

1/2 



(Ax) 



hJp 



■n 



1 + 1/1 + 



n + 2 



+ 2 



a . 



2 J Vyi + n/2 

which in the case of three spatial dimensions gives (Ax) m i n = 2.23533 hy/P . 

The states given in Eq. f)13[) . are normalizable states whose closure relation reads: 

^2 \ a 



(14) 
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(15) 



which turns out be a finite function for finite /3. The normalization constant iV can be chosen so 
that 5 n (£ — ^') reduces to the Dirac-5 distribution 5 n (£ — £') in the limit of zero minimal length 
(P — * 0) i.e. N = (2nn)~ n l 2 . Explicitly one has (for three spatial dimensions, n = 3): 

-(3/2+a) 



<5 3 (€-0 = r(i + a) 



2 Q (v^l^-^l 

7T 3 / 2 (^) 3 V V£ 



J. 



3/2+a 



(16) 



In Fig. [T] the closure function 63 is plotted against its dimensionless argument. It can be easily 
verified that it satisfies the relation J <i 3 £c>3(£ — £') = 1. 

For the purpose of deriving the first order correction in the minimal length of the Casimir-Polder 
forces it will prove useful to expand the closure function, c.f. Eq. (I15p . in powers of P and in 
terms of the Dirac-<5 distribution. In order to properly do this one must go to the momentum 
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representation (see appendix A for details) and one finds for the present model and within the 
KMM approach: 



5Jr 



l + (^) 2 f V= + ...] 5 n {r-r'). 



(17) 



where k = 2 + a. We take the previous equation as the definition of the numerical constant k 
whose value will in general depend on the model considered 2 . 

It should be noted that the KMM procedure to construct maximal localization states works prop- 
erly only for the generalized commutation relations of Eq. ([1]). Detournay, Gabriel and Spindel 
(DGS) 18)] have proposed a better definition of maximally localized states, which is based on a 
minimization procedure on the subset of all physical states, and not just on the subset of squeezed 
states, and turns out to be suitable for more general commutation relations than those considered 
in Eq. (p}. It should be realized that adopting this method (Ref. 18]) would only slightly affect 



our conclusions, namely Eq. (|17p . The maximally localized states are in this approach given by 
Bessel functions. In the end result, c.f. Eq. (I17p . we can anticipate that expanding the closure 
function in powers of the minimal length (hyf/3) it would merely amount to a change in the value 
of the numerical constant k which can however be expected to be always of order unity. 

Let us discuss in further detail this point. In particular we would like to show that adopting the 
DGS procedure, see Ref. is] ], the resulting regular closure function 63 (r — r') although represented 
by a different analytic expression is numerically very close to that computed in the KMM procedure. 
Furthermore its expansion in powers of (3 is as given in Eq. (I17p but with a slightly lower value 



of the numerical constant n. In Ref. 



181 ] the authors show that the maximally localized states 



corresponding to the particular model being discussed here (c.f. Eq. (jlip ) are given (in the spectral 
representation) in terms of Bessel functions: 



Vv ML (*) 



C n — exp 



n 



(18) 



where z = \z\; v = ^ — 1, /j, = iryfjj/^f 7 ! and the normalization constant C n depends again on 
the number of dimensions that we are considering. These states are characterized by a minimal 
uncertainty which is given by: 



(Ax) 



3u,l 

V2 



(19) 



2 We emphasize that k is a numerical constant and not a free parameter. Within the KMM procedure for obtaining 
the maximally localized states its numerical value (within the model considered here) is fixed by the number of 
spatial dimensions only. 
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FIG. 1: Plot of the regular closure function (hy/]3) 3 #3 (£ — £') along an arbitrary direction of the vector £ — £' 
(both positive and negative), versus the dimension-less quantity x = v2|£ — {th/J$)'i ( a ) the dashed line is 



the result obtained (Eq. [TSJ) within the KMM procedure of Ref. [l7| , with a = 1 + y/l + 3/2 corresponding 
to three spatial dimensions (n = 3); (b) the solid line is the result [Eq. [24] of this work] of the calculation 
performed with the (DGS) approach of Ref. [18]. This plot should not mislead the reader. Both functions 
satisfy the integral relation J d 3 £' S^(^ — £') = 1. The importance of numerical difference between the two 
functions, apparently inconsistent with this constraint, is reduced, when doing the integral, by the x 2 factor 
from the differential d 3 (£ — £'), and indeed both closure functions do integrate to 1, as has been checked 
both analytically and numerically. 



j v i being the first zero of the Bessel function J v (x). For n = 3 we have v = 1/2 and j v \ = tt so 
that the minimal length derived in the DGS procedure is (Ax) m i n = 2.22144 h\f]3 which is indeed 
smaller than the minimal length derived in the KMM approach 3 . With these states the closure 
function becomes: 



d^z [M t z)]2 exp 



>2;/ 



(20) 



Using the well known expansion for the Bessel functions: 



3 Recall that the minimization in the DGS procedure is carried through within the set of all physical states while in 
the KMM approach it is restricted to the squeezed states. 
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TABLE I: Within the rotationally invariant model (Model I) of Eq. (JTTJ (see Ref.jla, 117]) we compare the 
numerical constant k and the minimal uncertainty as discussed in the text for both the KMM and DGS 



procedure of obtaining the maximally localized states. Recall that the constant a = 1 + + n/2 appears in 
the maximally localized states of the KMM procedure. We also give the same details for the direct product 
model (Model II) used in Ref. [3], which is reviewed for the reader's benefit in Appendix A. 



k (n = 3) 



(Ax) r 



(n) 



(Ax), 



(n = 3) 



Model I (KMM) 

Model I (DGS) 2 

n 

Model II 10/3 



2 + a 

On/2-1, l) 2 



2 + 1 + w 4.58114 

7T 2 

2 + — w 5.28987 
10/3 



4 

Jn/2-1,1 



■a« 2.2353 



— « 2.22144 
\/2 



r(i/ + i) V2 



// + i 



+ 



(21) 



the constant C n can be chosen so that in the limit of a vanishing minimal length (5 — * (or /i — » 0) 
the regular closure function 5 n (£ — £') reduces to the Dirac-<5 distribution 5 n (£, — ^'): 



r(n/2) 



n/2-l 



(22) 



(2vrfi) n / 2 VM, 

In addition it is possible to derive an expansion of the closure function of the type in Eq. [T7] also 
within the DGS procedure. One derives the following expansion in powers of j3 (see Appendix A 
for details): 

2 + (>,i)7™ 



<^IV#f>=«n«-0 



1 + 



(^) 2 V| + .. 



(23) 



thereby defining, also in the DGS procedure, a numerical constant k, namely: k = 2 + (j^i) 2 : /n. 

On the other hand by analytically carrying through the integration in Eq. (|20p one finds the 
closed expression (when n = 3): 

1 2Si(x) - Si(x - 2vr) - Si(x + 2vr) 



(24) 



4vr 4 (/ly^) 3 x 

where again x = v 7 ^ - $,'\/{hy/P) and Si(x) = dt sm(t)/t is the sine-integral function. One 
again can verify that 63 (£ — £') satisfies the relation J d 3 £53(£ — £') = 1. Figure CD shows a 
comparison of the closure functions derived in the KMM and DGS procedure. Table I summarizes 
some results within the two approaches for the model of the present work (model I). 

The above discussion shows in particular that even within the same model adopting the KMM 
procedure or the more appropriate DGS method (where the minimization of the uncertainty is 
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done over all physical states) will produce slightly different values of the numerical constant k, 
thereby justifying the use of the widespread KMM approach at least within this particular model. 



III. QED IN THE TEMPORAL GAUGE 



In order to approach the quantization procedure in presence of a minimal length we review 
briefly the derivation of the temporal gauge [A (r, t) = ap(r,t) = 0] propagator in standard QED. 
The defining equation for the photon propagator in QED is: 



DV(r,t;r',t>) = % - (0\T [A*(r, t) A*(r f , t')] |0>, 



(25) 



from which, applying through the classical equation of motion for the fields Ai(r,t) (Maxwell 
equations) in the temporal gauge, it is easily verified that: 

i d 2 



c 2 dt 2 



V 2 ) 5« m - V*V r 



D mk (r : t;r',t' 



he 



5(t-t')(0\ E e (r,t),A k (r',t') |0) , (26) 



where the electric field operator E(r,t) = —(l/c)dA(r,t)/dt has been introduced. Canonical 
quantization prescribes then the equal-time commutation relations [xsj ] : 



E l (r,t),A k (r',t) = 4tt ihc 5 lk 5 s (r - r'), 
and the differential equation for the QED photon propagator becomes: 

i d 2 



V 2 I s* m - v*v m 



D mk (r -r'-t- t') = 4Tr5 lk 5(t - t')S s (r - r') 



c 2 dt 2 

The solution is found by going into a Fourier representation 

D rj {k-uj)= [ dtd 3 re iu)t e- ik - r D rj (r,t), 



and is given by the well known result: 

D rj (k,oo) = - 



4tt 



5 rj 



(27) 



(28) 



(29) 



(30) 



(u/c) 2 - k 2 + i0 f (oj/c)< 
As regards the calculation of the effect of a minimal length in the Casimir-Polder intermolecular 
interactions, following [2^], it will be convenient to use the propagator in the mixed representation: 

"+00 



dte iut D l J(t,r), 



which is found upon integrating in fc-space Eq. (|30|) : 

d 3 k 



kr 



(27T) J 



D ij (k,t 



(31) 



(32) 
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It is easily verified that D l3 (uj,r) satisties the following differential equation: 

,2 



which is solved by [201 ]: 



JS + V 2 ) 5^ - vV m 



D mk (uj,r) = An8 tk 5 A (r), 



Ik X 3i 



D*3(u,r) = 1 + 



IC c 



Mr uj 2 r 2 



+ 



x l x k ( 3c 2 3ic 



(33) 



r 2 I w 2 r 2 Mr 



,i\u\r/c 



(34) 



Quantization in minimal length theories 

Let us now discuss the quantization procedure of the electro-magnetic field in a quantum world 
with a minimal length. We shall proceed following the scheme adopted in 16|. In this case 



the procedure of canonical quantization gets modified because it turns out that the equal-time 
commutation relations of the fields, see Eq. (|27h . are different. Indeed now instead of expanding 
the field operators over plane waves (position representation wave functions of momentum states) 
one is forced to expand the fields over a set of maximally localized states \ip^ L ) of average position 
r (quasi-position representation of momentum states): 



a(p,X)e l (p,X) <^ /Z »e-^ 

+at(p,A)(^(p,A))*(p|^ /L )e + M , (35) 



introducing creation and annihilation operators which do satisfy the usual commutation relations, 



[a(p,A),aV,A')] = (2nf 5 X ' X ' 5 3 (p-p'), 



(36) 



with all other commutators vanishing (recall that momentum coordinates are commuting in the 
model which is being discussed here). The equal time commutation relations for the fields, c.f. 
Eq. P7|) . are then easily found to be modified to: 



E E (r,t),A k (r',t) = 4tt %hcS Us (VC |^ > 
= 4irihc5 ek 5 3 (r - r') , 



(37) 



so that, comparing with Eq. (|33p . the photon propagator in the presence of a minimal length, 
D' l:) (LO,r), is found from: 

,2 



D mk (uj,r) =4n5™6 3 (r) 



(38) 
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This approach to field quantization bears some similarity to that discussed in Ref. [2l[ where 
within a model of non-commutative space-time the real scalar field is defined as a mean value over 
coherent states. The non-commutativity is reflected by a modification of the Fourier transform 
replacing ordinary plane waves by gaussian wave-packets. We remark that in Ref.[21[ the authors 



derive for the Green function a differential equation which is similar to our Eq. ([38]) where the 
three-dimensional Dirac 5-function is replaced by a gaussian dumped by the non-commutative 
parameter which plays the role of the parameter /3 of this work. 

Being 5s(r) a regular function, it turns out that the standard QED photon propagator D r j(u,r) 
(from Eq. (|33j) ) is the Green function of Eq. (|38l) . the differential equation describing the photon 
propagator in minimal length QED, D l;1 (uj,r). Therefore it follows the central result of this work. 
In minimal length QED the photon propagator in the mixed representation is given by a convolution 
of the regular closure function 5s(r) with the standard QED photon propagator: 

D ik (u,r) = J d 3 r'D ik (uj,r - r')S 3 (r') . (39) 

From Eq.ll7lone can obtain, by means of partial integration, an expansion of the photon propagator 
in terms of the minimal length hy/J3: 

D ik (u, r) « D ik (u, r) + (H-\/]3) 2 ^ V 2 D lk {uo, r) . (40) 

Using the explicit form of the QED temporal gauge propagator given in Eq. ([34")) it is easily verified 
that: 



V 2 D lk (io, r) = -^-D ik {uj, r) , (41) 
c z 



so that we finally get: 



D ik (uj,r) 



1 " ? (^) 2 4 



D ik (io,r). (42) 



2 v V r-y c2 

A remark is in order at this point. It should be clear to the reader that the above approach to field 
quantization within a minimal length model breaks Lorentz invariance. This fact appears as well 
in other extensions of the standard model (SM) such as non-commutative quantum field theory, 
being al way s one of the major sources of debates between different authors. It should be noted that 
in Ref. 2lf] the question of Lorentz covariance is discussed and, within their non-commutative 2D 



model, successfully addressed. While completing this study the author became aware of a recent 
work [22j where a relativistic generalization of the Kempf algebra is proposed. Probably, based on 
this new deformed algebra, it will be possible to define a Lorentz covariant field quantization which 
would account for a fundamental minimal length. Certainly this point deserves further study, but 
goes beyond the scope of the present work. 
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IV. CASIMIR-POLDER INTERMOLECULAR INTERACTIONS 



The interaction of two neutral atoms (molecules) at rest in r\ and r 2 and with electric dipole 
moments d\ and d 2 is described by the operator: 



V = -E(r u t) ■ di(t) - E(r 2 ,t) ■ d 2 (t) 



(43) 



where E(r, t) is the operator describing the electric field. It is well known that, employing standard 
perturbation methods of quantum field theory, the interaction of the neutral atoms (or molecules) 
at distances r = |n — r 2 \, which are large compared to the atomic and/or molecular dimensions 
a: r 3> a, 4 is described in terms of (i) the dynamic polarizability tensor ot lk (uj) of the atoms: 



dr{0\T d\r), d k (0) |0)e 



+IU1T 



(44) 



n 



and (ii) the photon propagator in the mixed representation D tk (u>,r) [20(]. For the photon propa- 
gator we will use in the following the dyadic (tensor) notation D(u;,r). Assuming that the polar 
izability tensors of the two atoms are isotropic a$ 2 )( w ) = " a (i,2)i UJ )^ one obtains: 



OL\UJ 



+ 



3 ^ ' \iO n o — CO — i0 + LO n o + UJ — i0 + 

and the final expression for the interaction potential can be cast as: 

ih f +0 ° duj 



U{r) 



■uj ai{u)a.2{ijj)'Yr D(cj, r) ■ D(co, r 



(45) 



(46) 



2^7-00 2^ 

When a minimal length is present in the above expression one should replace the photon propagator 
D(u),r) with the modified propagator D(w,r). Using the approximate expression in Eq. (I42jl one 
finds: 

,21 



Tr 



D(w,r) • D(u,r) 



Tr 



D(u,r) ■ D(w,r) 



(47) 



so that the interaction potential is given by the usual standard model result with an additional 
term describing the effects of the minimal length: 

"+°° duj 



SU(r) = - % -^{hJp) 2 K J ^u;"n : uin 2 UMTr DU.r) ■ DU.r) 



By using the explicit expression in Eq. (|34|) it is easily found : 



Tr 



D(u,r) ■ D(w,r) 



1 + 



2ic 



5(f 



6ic d 



\uj\r co 2 r 2 



+ 



3c 4 



M 3 r 3 ' cu 4 r 4 



a 2i\uj\r/c 



(48) 



(49) 



4 In this field theoretic approach the distance r can otherwise be either much smaller (short distances) than the 
characteristic wavelength Ao, of the spectra of the interacting atoms (or molecules), comparable or much larger 
than Ao, (long distances). 
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The general expressions in Eq. (I46p and Eq. (I48h can be simplified in the limiting cases of short 
distances (a -C r <C Ao) and long distances (r 3> Ao). 
Short distances 

When r <C Ao, in the integral of Eq. (|48jl are important the values of w ~ loq ~ c/Ao so that 
cor/c -C 1 and in Eq. (|49[) one is allowed to keep only the last term and to approximate the 
exponential with 1. Performing the straightforward computations one then finds an interaction: 



(50) 



which scales as r~ 6 and is composed of two terms: (i) the well known London type potential |2C 
and the correction term which is due to the minimal length. In this regime of short distances the 
minimal length correction affects only the strength of the interaction and does not change its power 
law. 

Long distances 

In the limit of long distances r > Ao in the integral of Eq. (|48p are important only the values 
of a; < c/r <C luq. When u > u>q, the strongly oscillating complex exponential will suppress the 
integral. It is then possible to substitute the dynamic polarizabilities ct\ 2(0;) with their static 
values ai(0) and 02(0). We then find a correction of order (/ly^/r) 2 to the well known r -7 
Casimir-Polder interaction: 

Jg(r) = -j»«fMV»gsWaW, ( 51 ) 

8tt \ r J r' 

Thus, in theories with a minimal length, a correction term arises in the Casimir-Polder atomic and 
molecular interactions at large distances. The correction has the same sign {attractive potential) 
but scales with a different power law (r -9 ) relative to the QED result: 



, x _23 hca 1 (0)a 2 (0) 

^Casimir— Polder y" ) ~ ■ 



129 (WIX 
56 \ r J 



(52) 



4tt r 7 

This might be compared with the result of Ref. [l(| c.f. Eq. ([2]) which describes the Casimir energy 
of the plate-plate system. One can notice that while in the plate-plate case the minimal length 
correction derived in Ref. [lf| is repulsive and opposite to the standard result, in the case treated 
here of the interaction of two neutral atoms or molecules the minimal length correction is attractive 
i.e. of the same sign of the standard Casimir-Polder result. We should however be very careful 
in comparing the calculation of the Casimir effect of Ref. [r| with the present one. Indeed in 
Ref. [16(] a different model of generalized uncertainty relations is explicitly taken up. As opposed 
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f. Q g i 



to our Eq. (fTT|) the model proposed in Ref. [la] given by f(p 2 ) = l+Jh> 2 and g(p 2 ) = is the naive 



generalization of the one-dimensional model. Such model is known [23J] to be inconsistent with the 
KMM construction of maximally localized states. Due to the resulting non-commutativity of the 
coordinates the resulting n differential equations (corresponding to Eq. (|7|)) of the squeezed states 
cannot be solved simultaneously [23]. A different way of extending the one-dimensional GUP of 
Eq. ((I]) is to take a direct product of it. This so-called direct product model (model II) clearly breaks 
rotational invariance as opposed to the model of this paper (model I). The maximally localized 
states are obtained by taking the product of the one-dimensional states, see Eq. ([8]), along the 
different dimensions. This model is therefore characterized by different maximally localized states, 
which will in turn produce different closure functions which, when expanded in powers of f3, will 
provide different values for the numerical constant k as defined by Eq. (|17p (though we expect it 
to be always of order unity). In the appendix we provide details of the closure function within this 
model and deduce the corresponding value of k which turns out to be k = 10/3. We may therefore 
conclude that the Casimir-Polder minimal length correction within the direct product model is also 
attractive. That the Casimir effect for parallel plates within the direct product model (model II) 
as calculated in Ref. [lf| turns out to be repulsive may be due to the fact that the closure function 
there appears to have been computed not with the maximally localized states, but with the formal 
position eigenstates which are not physical states. Setting this point definitively calls for a detailed 
analysis of the Casimir effect using the proper maximally localized states within either model I or 
II and goes beyond the scope of the present work. 

The so-called Casimir-Polder force, for the atom-atom configuration, was first derived in 
1948 [24], but it was not measured definitively until 1993 [25], by looking at the deflection of 
an atomic beam passing through two parallel plates. For earlier measurements related to the the 



atom-plane configuration, see [2fJ|. There has been then a consistent renewal of interest with an 
increase of the measurement's precision as in Ref. [27] where this tiny interaction between a neutral 
system and a surface is reported. The authors employed a cloud of ultra-cold atoms in a Bose- 
Einstein Condensate (BEC) state. The range of distances explored was between 6 and 10 microns. 
See the recent review {28] for more details about the first generation of experiments, both in the 
Casimir effect and the related Casimir-Polder interactions, and for later developments and future 
directions. 

It seems that a better direction to explore the Casimir-Polder interactions experimentally is 
to consider the plate-s phe re (plate-atom) configuration. All the experimental advances have been 



achieved for this case 



281 ] . This suggests the direction where the present theoretical work could 
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develop: the calculation of the Casimir-Polder (plate-atom) interaction corrective term within a 
minimal length theory. 

V. DISCUSSION AND CONCLUSIONS 

We have studied the implications of models based on generalized commutation relations i.e. 
with a minimal length in the Casimir-Polder intermolecular interactions. The calculation is done 
following standard perturbation theoretical methods of quantum field theory and in particular the 
approach of I. E. Dzjalosinskij as illustrated in quantum field theory textbooks [2(j] is used. Here the 
interaction energy of two neutral atoms and/or molecules is related to the dynamic polarizability 
tensors of the two neutral bodies and to the photon Green function (propagator). The computation 
of the correction term due to the minimal length is thus carried out by discussing the QED photon 
propagator in presence of a minimal length. Quantization of the electro-magnetic field in the 
temporal gauge is performed in analogy to the canonical quantization procedure, the essential 
point being that the field operators A l (r,t) instead of being decomposed over a complete set of 
plane waves (momentum eigenfunctions) are now decomposed over a complete set of maximally 
localized states. This approach has also been followed by the author of Ref. [16| in deriving the 
Casimir potential energy of the plate-plate system in the presence of a minimal length c.f. Eq. ([2]). 
The main point is that the equal-time commutation relations of the field operators instead of being 
given in terms of a Dirac 5-function, c.f. Eq. (|27p . are now expressed by a finite regular function 
c.f. Eq. (I37p . and the photon propagator in minimal length QED is given by a convolution of the 
standard QED propagator with the regular function S^(r — r'). In order to compute the lowest 
order correction term the minimal length QED propagator is related to that of standard QED by 
performing an expansion in the minimal length parameter h^/J5. 

We have derived a corrective term to the (long distance) retarded Casimir-Polder interaction of 
two neutral atoms separated by a distance r 3> Ao, finding, c.f. Eq. (|52p . a new interaction term 
whose potential scales like r -9 as opposed to the standard r~ 7 result. 

Clearly should the minimal length (Ax) m i n = h\f]3 be of the order of the Planck length Lp 



the observability of this effect would be out of question. In Ref. [29( the study of the harmonic 
oscillator and the Hydrogen atom allowed to derive an upper bound on the minimal length by 
comparing with precision measurements on hydrogenic atoms and for electrons trapped in strong 
magnetic fields. The upper bound obtained is: 

(Ax) min = hy/p < 10-Hm. (53) 



16 



Assuming that the constant (5 appearing in the deformed Heisenberg algebra is a universal constant, 
this upper bound would presumably preclude a possible observation with Casimir- Polder interaction 
measurements with accessible distances r which are typically in the range 80 nm < r < 10 [im 



28 ]. Thus we conclude that taking into account the upper bound of Ref. [29|], the term derived in 



this work, c.f. Eq. (I52|) . would provide a correction whose relative strength we expect to be in the 
range: 

21 



10 



-v< s JL = o 
~ u 



(¥)' 



< io 



-18 



(54) 



which is clearly beyond any foreseeable improvement in the precision of measurements of Casimir- 
Polder interactions. As discussed in Ref. [30(] the rather strong upper bound in Eq. ([53]) could be 
avoided by assuming that the parameter (5 is not a universal constant and could vary from a system 
to another depending for example on the energy content of the system (the mass of the particle for 



instance) or the strength of some interaction. Indeed in Ref. 



301 ] by making this hypothesis the 



authors, through a comparison with the experimental resu 
in a gravitational quantum well (GRANIT experiment) 31 



ts for ultra-cold neutron energy levels 



32, 



33J derive a relaxed upper bound 
to the minimal length which turns out to be of the order of a few nano-meters [(Ax) m in < 2.41 
nm] . Even if this situation should apply, the relative strength of the corrective term will fall in the 
range: 

21 



io-<f = o 



< io 



-4 



(55) 



which still is out of the reach of current experiments in the Casimir-Polder interactions, whose 



precision is typically of 10 or 1% 



28j . With an increased experimental sensitivity to smaller 



length scales, of the order of a few nano-meters, and better experimental accuracy possible perhaps 
in the near future, the calculations described in the present work could find interesting applications 
and would provide a valuable tool to study the intimate structure of space-time. 
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APPENDIX A: DETAILS OF THE CLOSURE FUNCTIONS 



In this appendix we give details of the series expansion in powers of (5 of the closure function 
both in the model taken up in this work and in the one considered in Ref. 16] . In particular we 



show that in all considered cases the following relation holds with different values of k: 



s n (t -0= + | (^) 2 v? + 0(/? 2 )j 6 n (£ - £ 



(Al) 



It turns out that in order to get the correct series expansion in powers of j3 of the closure function 
one must resort to the unbounded momentum representation, while both Eq. 1151 and Eq. 1201 based 
on the spectral representation (compact space), can only reproduce correctly the first term of the 
expansion corresponding to j3 = 0. Let us consider this in detail within: (i) the rotational invariant 
model object of this work, (model I) detailed both within the KMM and DGS procedure and (ii) 
the direct product model of Ref. [16], (model II). 

Rotational Invariant Model (Model I) 

In Eq. (|15p we first perform a change of variables from the spectral representation to the momentum 
representation: 

v/1 + 2(ip 2 - 1 



with 



P 



(A2) 



The Jacobian matrix is: 



dzi 

— = 5ija + 2piPjb, 



where: 



y/1 + 2(3p 2 - 1 



PP 2 ' Pp a ^1 + 2(3p 2 

The Jacobian of the transformation J turns out to be given by: 

dzi 



J = det 



d Pj 



a n - 1 (a + 2bp 2 ) 



y/1 + 2(5p 2 - 1 

p P 2 



a) KMM procedure for the maximally localized states 



Eq. (rr5|) becomes: 



d n p ( a/1 + 2/3p 2 - 1 



n+a 



{2-Khy 



v/T+W 



exp 



i (c - a/1 + 2(3p 2 - 1 
ft«-0-P (Ip 2 



(A3) 
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We remark that the integration is unbounded. The integrand can be expanded in a power series 
around (3 = 0: 



(2vr/i)' 



i--M !Lt | ±£ p 2 -^«-0-pp a )+o(^ 



exp 



(A4) 



In the above expression one can make the substitution p — > iKV^ and thereby perform the mo- 
mentum integration obtaining: 

2 + a 



l-^-(^) 2 Vl + 0(/3 2 ) 



(A5) 



thus defining for this case the value of the numerical constant k to be given by k = 2 + a. We note 
that the last expression has been obtained after using identities (due to partial integration) valid 
for distributions such as: 



[(* " O • V € ] - a = -n V|<P(£ - o 



(A6) 



b) DGS procedure for the maximally localized states 

Using the maximally localized states of the DGS procedure one finds in the momentum represen- 
tation: 

1 [Ju{^z)f 



<U£-0 = c 2 n / d n P 



2 I ~ I Vl+W-1 



[3p 2 



y/1 + 2/V z 



2u 



Z -JZ-Z 



. Vl + 2/3p 2 -l 



exp 



i , yi + 2/y _ i 



(A7) 



where again the momentum integration is unbounded and the integrand admits a well defined 
series expansion in powers of (3. Using Eqj21] one finds similarly to Eq. | 



(2vrfr) r 



,. | n + 2 + fe, 1 )Vn p2 _| K _^ ) . pp2|+()| ,, 



exp 



- h ((-(') P 



(A8) 



and performing the momentum integration as before we end up with: 

2+(j^) 2 /n 



1 + 



(nv^) 2 v| + o(/3 2 



<5"(£-£'), 



(A9) 



therefore defining for the DGS procedure a value of the numerical constant k given by k = 2 + 
Uu,i) 2 /n. 
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Direct Product Model (Model I ) 

We note that while in Ref. [l6] the author announces to study the model: 



[xuPj] = ih(l + (3p 2 )5ij 



l,...,n, p 2 = p-p = ^2(p. 



2 



(A10) 



then he uses, as maximally localized states, the direct product of the one-dimensional states. 
More precisely the proper definition of the direct product model (model II) is {23)]: 



[xi,pj] = ih(l + 0Pi)Sij 



In this model the closure function is given by: 



d n p 



n 



cxp 



which factorizes in the product of n identical integrals: 



i = 1, . . . , n. 



J it c m Arctan(V/%) 



(All) 



8n(t-o = n 



n 



dp 



, 2ttH (1 + /3 P 2)2 
1 1 sin [xj] 



exp 



(A12) 



with Xj = vr(^- - £j)/(2fi>//J), while in Ref. Q it is found: 

'Trfe-^r 



(A13) 
(A14) 



- o = n 



sm 



n 



1 1 sin[xj] 



Wo remark I hat in Ref. 1.0]. in addition to some notational inconsistencies, the closure function 
appears to have been calculated improperly, as the (1 + (3p 2 ) factor appears in the denominator 
only with one power, while it should be squared, as one such factor is due to the change in the 
momentum measure and the other is from the square of the wave functions of the maximally 
localized states, see Ref. [15(. By performing a series expansion in powers of /3 of the integrand in 
Eq. (|A12p we end up with the following: 



l + |(^) 2 Vl + 0(/3 2 



(A15) 



which defines for this model n = 10/3, as reported in Table I. 

Let us remark that within this model, given by the direct product of the one-dimensional algebra 
discussed in section II, the KMM procedure gives the exact maximally localized states and so there 
is no need to distinguish between the KMM and DGS procedure [la ] . 
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It was however important to check that the sign of the constant k is the same as in model 
I. Thus even within the direct product model the Casimir-Polder correction due to the minimal 
length turns out to be attractive. 



[1] F. Wilczek, MIT internal note. 

[2] H. B. G. Casimir, Proc. Kon. Nederland Akad. Wetensch. B51, 793 (1948). 

[3] R. L. Jaffe, Phys. Rev. D 72, 021301(R), (2005); P.W. Milonni, Phys. Rev. A 25, 1315 (1982); J. 

Schwinger, Lett. Mat. Phys. 1, 43 (1975); C. Cohen-Tannoudij , J. Phys. (Paris) 43, 1617 (1982). 
[4] S. K. Lamoreaux, Phys. Rev. Lett. 78 5 (1997). 

[5] G. Bressi, G. Carugno, R. Onofrio, G. Ruoso. Phys. Rev. Lett. 88 (2002) 041804 

[6] C. Braggio, G. Bressi, G. Carugno, C. Del Noce, G. Galeazzi, A. Lombardi, A. Palmieri, G. Ruoso, 

D. Zancllo, Europhys. Lett. 70 (6), pages 754-760, 15 June 2005; W.-J. Kim, J. H. Browncll and R. 

Onofrio, EPL 78 21002 (2007); Woo-Joong Kim, James Hayden Brownell, and Roberto Onofrio, Phys. 

Rev. Lett. 96, 200402, (2006); V. V. Dodonov and A. V. Dodonov, J. Phys. A: Math. Gen. 39 6271-6281 

(2006); 

[7] E. Buks and M. L. Roukes, Phys. Rev. B 63 033402 (2001). 

[8] H. B. Chan, V. A. Aksyuk, R. N. Kleiman, D. J. Bishop and F. Capasso, Phys. Rev. Lett. 87 211801 
(2001) 

[9] J. C. Long, H. W. Chan, and J. C. Price, Nucl. Phys. B 539, 23, (1999). 
[10] A. Salam and J. Strathdee, International Centre for Theoretical Physics preprint IC/70/38 (1970). 
Based on a lecture given at the Symposium on "TTie Last Decade in Particle Theory" at the Center 
for Particle Theory, University of Texas at Austin, USA, April 14-17, 1970. Lett. Nuovo. Cimento, 4S1 
(101-180), 1970. 

[11] N. Arkani-Hamed, S. Dimopoulus, and G. R. Dvali, Phys. Lett. B 429, 263, (1998). 

[12] S. Hossenfelder, Calss. Quantum Grav. 23, (2006), 1815-1821; S. Hossenfelder, et al. Phys. Lett. B 

575, 85 (2003) [arXiv:hep-th/0305262| ; U. Harbach and S. Hossenfelder, Phys. Lett. B 632, 379 (2006) 

|arXiv:hep-th/0502i42] . 

[13] Kh. Nouicer and M. Debbabi, Phys. Lett. A 361 (2007) 305-315; K. Srinivasan, L. Sriramkumar and 
T. Padmanabhan, Physical Review D 58, 044009 (1998); Roberto Casadio, Alessandro Gruppuso, 
Benjamin Harms, Octavian Micu Phys. Rev. D 76, 025016 (2007). e-Print: larXiv: 0704 .22511 [hep-th] 

[14] P. K. Townscnd, Phys. Rev. D 15, 2795, (1977); D. Amati, M. Ciafaloni and G. Veneziano, Phys. Lett. 
B 216, 41 (1989); K. Konishi, G. Paffuti and P. Provero, Phys. Lett B 234, 276 (1990); M. Maggiore, 
Phys. Lett. B 304, 65-69, (1993); M. Maggiore, Phys. Lett. B 319, 83(1993); M. Maggiore Phys. Rev. 
D 49 , 5182 (1994); L. J. Garay, Int. J. Mod. Phys. A 10, 145 (1995). 

[15] A. Kempf, G. Mangano and R. B. Mann, Phys. Rev. D 52, 1108 (1995). 



21 



[16] Kh. Nouicer, J. Phys. A: Math. Gen. 38, (2005) 10027-10035. 

[17] A. Kcmpf, G. Mangano, Phys. Rev. D 55, 7909 (1997). 

[18] S. Dctournay, C. Gabriel, and Ph. Spindel, Phys. Rev. D 66, 125004 (2002). 

[19] C. Itzikson and J. B. Zubcr, Quantum Field Theory, McGraw-Hill, New York 1980 ISBN 0-07-032071-3, 
Chap. 3 pag. 127-128. 

[20] V. B. Berestetskii, E. M. Lifshitz and L. P. Pitaevskii, Quantum Electrodynamics, §85. Second Edition: 
Volume 4 of: "Course of Theoretical Physics" (Paperback) Butterworth-Heinemann, Oxford, 1982, 
ISBN: 7506 3371 9 

[21] A. Smailagic and E. Spallucci, J. Phys. A: Math. Gen. 36 (2003) L517-L521. 
A. Smailagic and E. Spallucci, J. Phys. A: Math. Gen. 37 (2004) 7169-7178. 

[22] C. Quesne and V. M. Tkachtuk, J. Phys. A: Math. Gen. 39 (2006) 10909-10922. 

[23] Toshihiro Matsuo and Shibusa Yuuichirou, Modern Physics Letters A 21 (2006), 1285-1296. 

[24] H. B. G Casimir and D. Polder, Phys. Rev. 73, 360 (1948). 

[25] C.I. Sukenik, M. G. Boshicr, D. Cho, V. Sandoghdar, E. A. Hinds, Phys. Rev. Lett. 70, 560 (1993). 
[26] Peter H. G. M. van Blokland and J. Theodoor G. Overbeek, J. Chem. Soc, Faraday Trans. 1, 1978, 
74, 2637 DOI: 10.1039/F19787402637 

The first experimental verification of the purely Casimir effect (in the plate-plate configuration) was 

reported in: M. J. Sparnay, Physica (Amsterdam) 24,751, (1958). 
[27] D. M. Harbor, J. M. Obrecht, J. M. McGuirk, and E. A. Cornell, Phys. Rev. A 72, 033610 (2005). 
[28] R. Onofrio, New Journal of Physics, 8, (2006), 237. 

[29] L. N. Chang, D. Minic, N. Okamura and T. Takeuchi, Phys. Rev. D 65, 125027, (2002) 
S. Benczik, L. N. Chang, D. Minic and T. Takeuchi, Phys. Rev. A 72, 012104 (2005); 
Earlier etimates relative to the Hydrogen atom had provided (Ax) m i n — h^fji < 10 -2 fm. 
See F. Brau, J. Phys. A 32, 7691, (1999). 

[30] F. Brau and F. Buisseret, Phys. Rev D 74, 036002, (2006). 

[31] V. V. Ncsvizhevsky ct al. Nature (London), 415, 297 (2002). 

[32] V. V. Ncsvizhevsky et al. Phys. Rev D. 67, 102002, (2003). 

[33] V. V. Ncsvizhevsky et al. Eur. Phys. J. C 40, 479, (2005). 



